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1. Introduction 

QCD is a strongly coupled theory, where nonperturbative physics plays a crucial role. As such, 
it is hard to handle analytically. A way out to investigate nonperturbative physics is model building 
(e. g. by using holographic QCD models, effective models like the PNJL model or others). Another 
extremely powerful tool is simulating QCD on a finite lattice. One can also attempt to quantize the 
theory in the continuum and try to get as good as possible information out of this by a variety of 
techniques. The latter philosophy shall also be employed in this work. 

With lattice QCD, expectation values of gauge invariant operators can be computed without the 
need of gauge fixing. The eventual numerical estimates for physical quantities like a particle's mass 
are in good agreement with experimental data. The merit of lattice QCD is that it can also provide 
us with physical information in theories which do not appear in nature. Some famous examples are 
QCD with a number of colors other than 3 or QCD without quarks (gluodynamics). In the latter 
case, the physical spectrum supposedly exists of colorless pure glue states, the glueballs, see ||] 
and references therein. 

We recall the classical SU (N) Yang-Mills action in d = 4 Euclidean space time, 

S YM = \j f 'd 4 xF« v F£ v . (1.1) 
This action possesses an enormous local invariance w.r.t. 

A M -> A s ^ = S + d^ l S + S + A fl S SeSU(N), (1.2) 

or in infinitesimal form 

A*+Df(O b , Df^d,8 ab - g f bc A c ^. (1.3) 

We need to reduce this enormous overcounting of physically equivalent gauge configurations by 
fixing a gauge. In principle we have a complete freedom to do so. Usually, one can pick a gauge 
suitable for the problem under study. However, there is an important restriction, as one should be 
assured that the eventual gauge fixed theory needs to be renormalizable. Not every gauge belongs 
to the class of renormalizable gauges. 

We shall now focus on one particular example of a renormalizable gauge, viz. the Landau 
gauge, d^A^ = 0. This is a very popular gauge in the continuum, as it has many nice (quantum) 
properties ph. According to the Faddeev-Popov procedure, the gauge fixed action reads 

Sym + S gf = J d 4 x fe v + b a d^ +l?d ll DffA . (1 .4) 

This gauge fixed action no longer enjoys a local gauge invariance. However it gets replaced by the 
equally powerful nilpotent BRST symmetry, s(Sym + Sgf) = 0, 

sAl = -Dfc\ sc" = ^f ahc c b c c , sc a = b a , sb a = 0, s 2 =0, (1.5) 

which can be used to prove the perturbative renormalizability and perturbative unitarity of the 
theory. 
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2. Gauge (Gribov) copies and the Gribov-Zwanziger (GZ) approach 

During the Faddeev-Popov procedure, it is always tacitly assumed that there is one and only 
one solution to the gauge fixing condition. Gribov was the first to realize this to be wrong and 
constructed explicit examples in his seminal work [ph. If we take in the Landau gauge, d^A^ = 
0, and consider an (infinitesimal) gauge transformation, A' =A^ +D^co, then we quickly see that 
dfiA'n = if d^D^co = 0. Apparently, we encounter gauge copies if the Faddeev-Popov operator 

M ab = -d ll D ab (2.1) 

has zero modes. In order to exclude these copies from the path integral, Gribov proposed to restrict 
the integration to the Gribov region Q. where dA = and M > 0. This CI corresponds to local 
minima of the functional / d 4 xAj 1 along the gauge orbits. This is already an improvement of the 
original Faddeev-Popov quantization procedure. The question turns out to be how to implement 
this kind of restriction to CI in the continuum formulation? Gribov and later on Zwanziger 
worked out this problem and proved many properties of the region CI, for example that every gauge 
orbit passes through CI We should however warn that CI still contains copies, not related to zero 
modes of the Faddeev-Popov operator [0]. 

After a lengthy analysis, Zwanziger was able to implement the restriction to the Gribov region 
CI to all orders by means of a local action, known as the Gribov-Zwanziger (GZ) action. The 
eventual GZ partition function becomes ^ 



with 



l FP = J [dA] [dc] [dc] [db] [dp] [dp] [do] [dw]fT s ° z , (2.2) 



Sgz — Sym + S g f + Sfpyaa + S r , 

S^am = f d 4 x(lp a ;d v (d v (p^ + gf bm A b v (p>» c ) -a%d v (d v (o^ + g f bm A b v co'; c 
-g(d v G>™)f bm (D v c) b (p™ 



Sy = -fg J d*x(f^(p* c +f bc Affi + ^ {N 2 -1)A , (2.3) 

where extra bosonic (p, p) and fermionic (ft), To) fields were introduced. The parameter / carries 
the dimension of mass and must be self-consistently fixed to the nonzero solution of the following 
gap equation 



vac 



dy -0, ,2.4) 

commonly known as the horizon condition [Sj, thereby giving y ~ Aqcd, a typical example of 
dimensional transmutation. It can be easily checked that for y = one recovers the Faddeev-Popov 
theory. It is important to mention that this action ( |2T3| ) defines a renormalizable theory see ^ |8|] 
and references therein. 

What about the BRST symmetry? One can naturally extend the BRST symmetry to the new 
fields 1 

*P£ C = < C , sco^ = 0, sW™ = p^ c , s7p a ' =0 . (2.5) 



For y = 0, we then obtain a trivial extension of the usual Faddeev-Popov gauge theory. 
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The symmetry of the original action under the BRST transformation is softly broken 

sS GZ =gfj d*x (f abc AlG>l c - (D«mc m )(lj% + <pj c )) . (2.6) 

With softly we mean that it is proportional to the mass parameter y 2 , thus it can be controlled at the 
quantum level. Very recently, an equivalent formulation of the GZ theory was given in which case 
the breaking is even converted into a linear breaking [g]. 

Apparently, treating gauge copies a la GZ leads to a loss of the BRST symmetry. As such, the 
situation of how to define physical states also becomes less clear (see later). 



3. Why studying propagators? 



During the past decade, a lot of effort went into the investigation of the elementary gluon 
and ghost propagators. One might wonder why so much study is devoted to these gauge variant 
quantities, as these do not correspond to physically measurable quantities? Propagators are the 
basic building blocks of quantum field theory: they do describe the propagation of the elementary, 
albeit perhaps unphysical, degrees of freedom and they are the "simplest" objects to compute. 
In any Feynman diagram-based approach to QCD, propagators appear. This reaches far beyond 
perturbation theory, one needs only to think about Schwinger-Dyson, Bethe-Salpeter, sum rules, 
. . . approaches to QCD, which are all nonperturbative in nature. In addition, since gluons are 
confined at low energy, we might expect to see something nontrivial already at the level of the 
propagators. 

These propagators were also intensively studied using lattice simulations, in particular in the 
Landau gauge. This gauge is very suitable to be simulated, as it corresponds to searching for the 
minima of the functional / d 4 .xA 2 along the gauge orbits. As such, one gets the "exact" propa- 
gators and can compare those with analytical results in order to test the latter. If analytical and 
lattice results are in good agreement, we may have a certain degree of confidence that the analyt- 
ical approximation scheme gives reasonably good propagators, which can then be used in other 
computations which depend heavily on knowledge of nonperturbative propagators. 



In Figure 1, taken from [ 10], we show the lattice SU(3) Landau gauge gluon propagator and it 



Renormalized Gluon Propagator - (I - 3 GeV 



a f 
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• (5 = 5.7 
■ (5 = 5.7 

J = 5.7 
- 3-5.7 
f> = 5.7 



Figure 1: Gluon propagator renormalized at jj. = 3 GeV. 
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is clear that there is no sign of a blow-up in the infrared, typically occurring when using perturbation 
theory in the massless Faddeev-Popov scheme. This is indicative of nonperturbative effects, and 
we shall try to motivate that these might be related to Gribov copies. 



4. The Refined Gribov-Zwanziger (RGZ) approach 



Using the action ( |2.3| ), it is readily verified that the tree level gluon propagator reads 

„2 



Dip 1 



P 4 + X 



4 ' 



(4.1) 



which is indeed infrared suppressed, although it vanishes, in contradistinction with the lattice Fig- 
ure 1. We have set here A 4 = 2g 2 N r / i . These effects persist upon including loop corrections [pT|]. 
For the ghost propagator G(p 2 ), one can prove at any order |Q] or compute explicitly [12] that by 
invoking the gap equation (E 



p 2 G(p 2 )~\, forp 2 ~0, 



(4.2) 



which again seems to be at odds with large volume lattice data [|13[]. 

Apparently, something is missing in the Gribov-Zwanziger formulation as standing. In order to 
overcome this, in [14, 15] extra dynamical effects due to nonperturbative d = 2 condensates were 
taken into account. We recall that the d = 2 condensate was popularized in the last decennium, 

1^, 18]. As a result of the analysis, one finds a gluon propagator of the 



thanks to works like 



form [14, 15] 



D( P 2 



p 2 +M 2 



p A + (m 2 +M 2 )p 2 + M 2 m 2 + 2g 2 N r f 



whereby m 2 and M 2 are mass scales corresponding to condensates, in particular 



m ~ 



{A 2 



M 2 ~ (qxp-coco) 



(4.3) 



(4.4) 



We observe that D(p 2 ) is still infrared suppressed, but D(0) ^ thanks to the presence of M 2 . 
Hence, the gluon propagator lattice data is already qualitatively reproduced. For the ghost, one 
finds [14, 15] G(p 2 ) ~ -y for p 2 ~ 0, again consistent with lattice data. 



Having found a reasonable qualitative agreement, one might wonder whether the lattice data 
for the gluon could also be quantitatively fitted with a propagator of the form (fO|). This was tested 
in JIfj|]. The fit with m 2 = did not work out well, indicating that (A 2 ) is of importance. The 
following continuum extrapolated values were reported 



M 



2.15±0.13 GeV 2 , 



1.81 ±0.14 GeV 2 , 



2g 2 Ny 4 



4. 16 ±0.38 GeV , (4.5) 



leading to the fits displayed in Figure 2 for the gluon propagator and its form factor. We see that 
the fit works out well up to p ~ 1.5GeV. As a byproduct of this analysis, an estimate for the d = 2 
gluon condensate was obtained, 



? 2 a 2 : 



ji=10GeV 



3GeV 2 , 



(4.6) 
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Renormalized Gluon Propagator - (i = 3 GeV Renonnalized Gluon Dressing Function - p. = 3 GeV 




Figure 2: Fit to the lattice gluon propagator (left) and form factor (right). 

which is in the same ballpark as other, rather independent, approaches, see [ |To| ] for references. 

It is apparent that the RGZ framework seems to be able to describe quite well the nonperturba- 
tive (gluon) propagator. One is thus lead to believe that RGZ can be a good starting point to study 
nonperturbative aspects of Landau gauge QCD. Future work should be devoted to try to compute 
the scales in D(p 2 ) in a clean analytical fashion using an effective potential approach. Recently, 
in JT9| ] an even more general RGZ setting was proposed. We conclude this first part of our talk 
by drawing attention to similarly good-looking results for the two basic propagators in other ap- 
proaches obtained by solving Schwinger-Dyson equations, see some of the other proceedings and 
PL Pi- 



5. Glueballs in the (R)GZ approach and the concept of /-particles 

There is however more to life than gluon and ghost propagators. As the latter are believed to 
be unphysical, one should look what the physical degrees of freedom could be, such that these can 
be described by the (R)GZ theory. Therefore, we need to look at glueball correlation functions. 
The big question is how to proceed from the "simple" gluon/ghost to the "less simple" glueball? 

It turns out to be useful to consider two sets of variables, being the old variables (gluons, 
ghosts and extra GZ fields), which are useful for renormalization issues, and a set of new variables, 



the so-called /-particles [22], useful for spectral issues 



Let us first introduce these /-particles. We notice that there are 2 complex (conjugate) gluon 



"masses" in the propagator (4.3), given certain bounds on mass scales, which are for example 
fulfilled by the numbers (4^5). A set of 2 cc masses immediately leads to a tree level gluon positivity 
violation (which is also seen on the lattice), meaning that the gluon cannot be a physical excitation 
[26]. One can believe this is reflective of gluon confinement [|7|]. We also observe that there is a 



field mixing in the tree level RGZ action 

/ d 4 x ( l A^(-d% v -d,d v )Al + r^-d 2 )^ b + fgf abc A^ . (5.1) 
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These cc masses are a bit hidden in these variables; we shall also set m 2 = for the remainder of 
this talk. By using a set of linear transformations the foregoing tree level action can however 
be recast in the following form 



d 4 x 




-d 2 + )^ + 2< 1-5 + 2 ] T ?M+ rest 



Clearly, the new fields A„ and 7]fJ have cc masses. They contain the gluon fields, which are un- 
avoidably mixed with the (R)GZ content due to y ^ 0. We call these novel fields the /-particles of 
the (R)GZ theory. The question how to describe glueballs in the (R)GZ context is still unanswered. 

A first potential pitfall is how to define a physical subspace without BRST invariance? These 
/-particles are obviously unphysical, which can be called "confinement" if one is of a very optimist 
nature. But gluon confinement is of course more. We should be able to define a physical subspace 

of purely gluonic states (i.e. glueballs), which should not decay into unphysical gluons/ghosts/ 

It looks like we need a symmetry to define such a subspace, thereby expelling the unphysical stuff. 
This smells like a BRST symmetry application, but the GZ action breaks BRST as we have seen 
already. In it was shown how to construct a local albeit not nilpotent BRST symmetry of an 
equivalent version of the GZ theory, but it remains unclear if this symmetry is sufficiently strong to 
define a physically sensible set of glueball operators. 

A second pitfall concerns the renormalization of a suitable glueball operator. As a glueball 
is a kind of bound state of gluons, we need a suitable (local) composite operator, whose quantum 
numbers correspond to the glueball state under investigation. As glueballs are physical, we expect 
gauge invariant composite operators 2 . Again, the loss of the quantum version of the BRST sym- 
metry seems problematic. Nevertheless, in case of the scalar glueball, it was shown in [ ]24| ] that a 
renormalization group invariant extension of the classically gauge invariant operator F 2 v (x) can be 
constructed in (R)GZ. As expected, there is operator-mixing into 

O = F 2 V + s(d = 4 operators) + J 2 (d = 2 operators) + EOM-terms . (5.2) 

Notice that, unlike in usual QCD, 



{0{x)0(y))^{F' v {x)F 2 a {y)) , (5.3) 

the reason being the BRST breaking ~ y 2 . This means that the mixing terms do influence the 
correlator. 

How should a general glueball operator 0(x) look like in the (R)GZ world? 0(x) should be 
renormalizable, and for y = 0, we expect to find back the original QCD cohomology output, thus 
something of the form 

Oy = o(x) = gauge invariant operator + s(d = 4 operator) + EOM-terms . (5.4) 

Hence, for y / 0, we seem to be driven to 

O r ^o(x) = gauge invariant operator + s(d = 4 operator) 

+Y 2 (d = 2 operator) + EOM-terms . (5.5) 

2 Or more precisely, quantum BRST cohomology classes, if we have a nilpotent BRST operator. 
3 EOM-terms stand for terms proportional to the equations of motion. 
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The discussion on renormalization of classically gauge invariant operators is most easily given in 
the old variables. The story changes when we want to check whether the operator generates a 
physical two-point function. In order to speak about a physical propagator A(p 2 ) = (0(p)0(—p)), 
D(p 2 ) must have a decent spectral representation 4 , 

A(p 2 ) = -^— 2 + f ^\dt , p 2 eC (5.6) 

so that A(p 2 ) displays a branch cut only along the negative axis, with positive discontinuity p(t), 
given by 

pW = t- lim [A(-t-ie)-A(-f+ie)] . 

2711 e^0+ 

This positivity can be easily understood from the optical theorem since p is proportional to the 
cross section, which ought to be positive. The m 2 > correspond to physical particle masses, 
while To > corresponds to the multiparticle-threshold. 



As demonstrated in Q22[] , the /-particles are very suited to derive the spectral representation of 
e.g. the correlation function built with F 2 V . At lowest order, we can work in the quadratic tree level 
(Abelian) approximation, in which case F 2 V , rewritten in /-particle field strengths, reads 

o(x) = jf%f{iv = ^v%v + 2^v^v + 2 ^vH^v • ( 5 - 7 ) 

A priori, it is unclear how to compute the spectral representation, as the Cutkosky cut rules are in 
principle only intended for use with real masses in Minkowski space. This problem can neverthe- 
less be handled, but discussing this here would lead us to far. Let it suffice to mention that a few 
examples were worked out in [fffi). 

For oi (x) = A^ V T]2 V , it turns out that 

( 0l (p)oi(-p))= f withp(0>0. (5.8) 

This is good news, as this represents a physical spectral representation. On the contrary, 02 (x) = 
1/2^ + 1/2^,1); leads to 



(o2(p)o 2 (-p)) = I tt4+cc (5.9) 

J curve in C I ~T p 

thereby displaying cuts in the complex plane. Unfortunately, this also means that the tree level 
version of F 2 V itself, viz. f 2 v = o\ +02, leads to an unphysical correlator [22]. 
It would then seem that only taking A^ V T]^ V is a good choice, since 

^v^v=/M 2 v+ reSt > (5-10) 



and it leads to a physical correlation function. However, this operator falls outside the class (p.5|), 
and it therefore seems to be doubtful to be renormalizable, i. e. controllable beyond the tree level. 
In addition, we do not have any information yet on the spectral properties beyond tree level either. 



We work in Euclidean space. 
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We are currently investigating an operator of the kind 

f\iv + s(other operators) + y 2 (other operators) , (5.11) 



restricted to the Abelian level. If the spectral properties would turn out to be OK, the form (5.11) 
would allow to at least write down an extension of the operator to the quantum level. One can 
then try to investigate its renormalization and, if possible, its higher order spectral properties. This 



kind of operator ( p. 1 1[ ) would also fit with the new BRST constructed in [|23|]. From the tree level 



results, if they are physically decent, one can also already extract information on glueball masses, 



perhaps along the lines of [ ]25[ ] where a GZ-like theory was studied in relationship with glueball-like 
operators. 

We conclude that it is apparently very hard to accommodate good renormalization and good 
analyticity properties at the same time when it comes to the study of glueball operators. We believe 
this is not only of relevance to work in the (R)GZ context, but to all people active in propagator 
research: how can one go from the unphysical gluon/ghost propagators to a well-defined physical 
subspace of glueball operators, which can be controlled at the quantum level. In addition, one 
should also try to find reasonable estimates for the glueball masses [Q]. This looks like a very 
ambitious program, but we hope that it will stimulate a lot of new research in the coming years. 
It would also be interesting to find out whether the difference between GZ and RGZ plays a role 
when it comes to glueball properties. Also the issue of a BRST symmetry needs to be further 
clarified, even in the Schwinger-Dyson context due to the potentially subtle role played by boundary 



conditions [20, E8L 29]. 
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